Chemistry in Evaporating Ices 



— Unexplored Territory 



Cesare Cecchi-Pestellini 

INAF - Osservatorio Astronomico di Cagliari, Strada n.54, Loc. Poggio dei Pini, 09012 

Capoterra ( CA ), Italy 

ccpOca. astro . it 

Jonathan M. C. Rawlings 

University College London, Department of Physics and Astronomy, Gower Street, London 

WC1E 6BT 

j crOstar . ucl .ac.uk 

Serena Viti 

University College London, Department of Physics and Astronomy, Gower Street, London 

WC1E 6BT 

sv@star . ucl .ac.uk 

David A. Williams 

University College London, Department of Physics and Astronomy, Gower Street, London 

WC1E 6BT 

dawOstar . ucl .ac.uk 
ABSTRACT 

We suggest that three-body chemistry may occur in warm high density gas 
evaporating in transient co-desorption events on interstellar ices. Using a highly 
idealised computational model we explore the chemical conversion from simple 
species of the ice to more complex species containing several heavy atoms, as 
a function of density and of adopted three-body rate coefficients. We predict 
that there is a wide range of densities and rate coefficients in which a significant 
chemical conversion may occur. We discuss the implications of this idea for the 
astrochemistry of hot cores. 



Subject headings: astrochemistry - ISM: molecules - ISM: clouds 
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Introduction 



Hot cores are well-known as sites of relatively complex chemistry in the interstellar 
medium 



I see, e.g., 



Snyderl 120061 for a discussion of Sgr B2 N-LMH). None of these species is 



readily formed in conventional gas-phase chemistry under typical conditions of interstellar 
clouds nor of those pertaining in hot cores (number densities nn 2 ~ 10 7 cm -3 , and temper- 
atures T ~ 200 — 300 K). Attention has therefore focused on interstellar ices as a potential 
source of these relatively complex species. Laboratory experiments have shown that irra- 
diation of ices of similar composition by fast particles or short wavelen gth electromagneti c 
radiation can induce greater chemical complexity to arise in the ices (e.g.. lOberg et al.ll2008l ). 
Concurrently, e xtremely detailed computational models by Herbst, Garrod, and their col- 



laborators (e.g., iGarrod. Widicus Weaver fc Herbstl 120081 ) predict that a rich and relatively 



complex chemistry can slowly arise in the low temperature ices. Then, when the ices are 
evaporated in the very dense and warm gas in the vicinity of a newly-forming star, these 
complex molecules are released to the gas phase and can be detected as hot core molecules. 

In this work, we suggest an alternative approach that may work in parallel with these 
chemistries to form the large organic molecules detected in hot cores. We ask: can three-body 
gas-phase chemistry in high density transient events on evaporating ices create complex 
molecules? In the process of becoming a hot core, material warms up over a finite period 
of time from a low temperature (around 10 K) to about 200 K; laboratory experiments 
( ICollings et al.ll2004l ) show that desorption occurs in several distinct and narrow temperature 
bands, of which the most important for our present purposes is the so-called co-desorption 
band. This band is when the major component of the ice, H 2 0, desorbs and carries with 
it all other species. This understa nding of desorpt ion has been shown to be consistent with 
current observations of hot cores ( IViti et al.l 120041 ) . 



In the picture we present here, we assume that during the warm-up of a pre-protostellar 
core a major part of the ices is abruptly converted during the co-desorption event from solid 
phase to gas phase. A truly instantaneous conversion from solid to gas would create a gas 
with a number density similar to that of the solid, i.e. about 10 23 cm -3 . However, it is 
unlikely that the gas attains such high densities, but it is nevertheless possible that the 
density is initially at a very high level indeed, although for a very short period of time. It 
is our purpose here to explore whether three-body reactions in extremely dense and fairly 
warm gas can create molecules with complexity similar to those observed to exist in hot 
cores. The source species available for these chemical syntheses are assumed to be those 
molecules available in unprocessed interstellar ices. Can we use H 2 0, CO, CH 4 , etc., to 
create molecules of greater complexity in a highly transient and extremely dense phase? 



Ideas of a similar kind were first explored by iDuleyl (120001 ) who suggested that amino 
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acids, peptides, and a variety of organometallic compounds could be created from evapo 
rating ices confined within cavities inside aggregate grains. Later ICecchi-Pestellini et al. 



( 120041 ) computed radiation field intensities within such cavities, and suggested that pre- 
biotic molecules might arise in the chemistry promoted in those cavities. The processes 
described in these papers would take place over long intervals of time during the lifetime 
of the core gas. In this paper, by contrast, we consider whether the conversion from chem- 
ical simplicity to chemical complexity can occur very rapidly within the transiently very 
high density gas arising during t he evaporatio n of ices in the co-desorpt i on ph ase. Of 
course, the processes discussed bv lDulev (120001 ) , by ICecchi-Pestellini et al.l (120041 ) , and by 
Garrod. Widicus Weaver fc Herbstl (120081 ). and the ideas presented here, are not mutually 
exclusive. 

In Section 2 we describe the physical and chemical models that we have adopted. The 
evaporating gas is assumed to be initially at a very high density (treated as a free parameter) 
and to expand either freely into a vacuum, or more slowly if there is some partial confinement 
of the gas by the dust grain geometry. The chemistry that we adopt is hypothetical, since 
there is little information available about three-body gas phase chemistry at exceptionally 
high number densities. Our work, therefore, is intended to be a feasibility study exploring a 
hitherto highly transient region of parameter space. We present in Section 3 the results of 
our computational modelling. Our conclusions are given in Section 4. 



2. The model 



2.1. The initial chemical conditions 



The co mposition of interstellar ices, measured along differen t sight-lines, varies signifi- 
cantly (e.g., Iwhittetll2002[ iGibb et al.ll2004[ iBoogert et al.l 120041 ) . We adopt a composition 
within the observed range. The precise composition is unimportant for our present purposes, 
where order-of-magnitude abundances of product molecules are adequate. We assume that 
the instantaneously sublimated ice mantle has a typical composition as given in Table [TJ 
For our feasibility study we currently ignore other species, whose abundances are less than 
1% of the mantle composition. We also do not include C0 2 , which has a typical relative 
abundance of 20%, on the assumption that it is tightly bound and unreactive in the assumed 
conditions. 
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2.2. The chemical network 

With this limited species set, the likely reaction channels are also somewhat restricted. 
As the composition of the evaporating gas is dominated by H 2 and, to a lesser extent, by 
CO we assume that the third body in the three-body reactions is H 2 0. 

Three-body reactions are believed to have important roles in astronomy. They deter- 
mine the chemistry in planetary atmospheres; for example, the reaction between molecules 

CH 4 + H 2 + M -> CO + 3H 2 + M 



(whe re M is a third body) may affect the deep water abundance on Jupiter (IVisscher et al. 
2010l ). The same authors also suggest that the reaction 



H 2 + H 2 CO + M -> CH 3 + OH + M 

may open the way to radical formation (in the present work, we have excluded radical 
formatio n and react ion). Similarly, three-body reactions play a role in cool stellar atmo- 



spheres ( jTsujil Il973l ) in the formation of simple species such as H 2 , CO, C0 2 , CH 4 , NH 
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scribed by 


Garrod. Widicus Weaver &; Herbst 


2008) 


1992: 


Tielens & Haeen 


1982; 


Allen & Robinson 


1977L 



These surface chemistries invoke rad- 
ical/molecule reactions in which the surface is the third body. It should be noted that all 
the above chemistries rely almost entirely on theoretical estimates rather than laboratory 
determinations. 

Our proposal is that three-body reactions between molecules (rather than radical/molecule 
reactions) may generate a product channel of astronomical interest. Such a channel would 
depend on bond-breaking and atomic rearrangement, and could be driven by the high energy 
tail of the Boltzmann distribution in the evaporating gas. Examples of this kind of reaction 
that might be possible are those quoted above for Jupiter, or an alternative channel 

CH 4 + H 2 + M ->• CH3OH + 2H + M 

possibly followed by 

CH3OH + CO + M -> CH3COOH + M 
An interesting pair of similar reactions is quoted in the NIST databasd^ 



H 2 + CO + M ->■ HCOOH + M 



1 http://webbook. nist.gov/ 
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followed by 

HCOOH + H 2 + M HOCH 2 OOH + M 

Thus, it appears that reactions of the type we consider here are under consideration in 
various applications. We therefore assume that three-body reactions may occur between 
molecular species such as those of the initial composition. 

We do not claim to know what the products of the various reactions are; we only 
hypothesise that non-defined complex organic molecules (which we label Pi to P26) may be 
formed. Obviously, this is highly simplistic and does not include the possibility of complex 
branching ratios, or the other (small species) products of the reactions. It therefore follows 
that we cannot track the abundances of the parent species with any degree of accuracy and, 
once the model indicates that their abundances are reduced through conversion into more 
complex molecules, then the model breaks down. 

We have considered various types of possible three-body reactions. Normally, one thinks 
of the process as one that involves a chemical reaction between two species. The third body, 
which is chemically inert, then collisionally stabilises the excited product of this reaction. 
However, we can also envisage reactions in which all three species are chemically active. The 
reactions in this category are listed in Table [2j As stated above, we assume that H 2 is 
always a partner in these reactions. 

The more usual formation channels would involve one of the three reactants as a chemi- 
cally passive partner. Thus to obtain the products Pi to P i5 above would require two stages 
of reaction, involving intermediate species (Pi6 to P 2 6) after the first stage. These, first stage, 
reactions are listed in Table [31 Again, the third (now passive) partner in these reactions is 
taken to be H 2 0. 

The second stage reactions, resulting in the formation of Pi to P 15 , are given in Ta- 
ble HI Note that we further assume that the sequence of reactions to form a product is not 
important; thus reaction 16 followed by reaction 27 is completely equivalent to reaction 22 
followed by reaction 47. 

This could be extended to include reactions between less abundant species (e.g., CH 4 
and OCN~ etc.) but the products of these reactions would obviously be less abundant than 
the products considered here. We shall ignore these minor routes. 

The rate coefficients for these various reactions are entirely hypothetical as no detailed 
information is available for any of the reactions in our list. We assume that the temperature 
dependences of the rate coefficients can be described by the usual Aarhenius formalism 



k = k (T/300 K)° exp (-0/T) 



(1) 
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For the purposes of this exploratory calculation we assume that all reactions have the same 
basic rate coefficients kn. By referring to existing databases of three-body reactions (e.g. 



RATE06, IWoodall et al.ll2007f ) it is evident that most of the values of ko for the reactions lie 
in the range 10~ 27 — 10~ 32 cm 6 s^ 1 . The values of a typically lie between and —3, whilst the 
j3 values are very variable, and are quite often negative for the temperature range for which 
a fit is given. We consider two possible values for a: and —2. So as to avoid excessive 
degeneracy in the parameters we set /3 to zero and vary k to incorporate the effects of a 
possible barrier. Any particular value of the rate coefficient can be re-interpreted in terms 
of a canonical value of k Q together with a particular value of (3 (see Section 4). 

These are the only reactions in our chemical network; at the densities that we are 
considering, photoreactions and cosmic-ray ionization reactions are entirely negligible and 
so have been omitted. 



2.3. The Physical Model 

We can estimate very roughly the expected timescale for significant conversion of re- 
actant species into complex organics on the assumptions that (i) there are no significant 
activation energy barriers to the reactions in question, and (ii) the reactions are generally 
constructive - that is to say at least one the products is more complex than the reactant 
species. If the density is taken to be n ~ 10 23 cm -3 , the rate coefficients are all 10~ 30 cm 6 s _1 , 
and the fractional abundance of the reactants (e.g., CH4, OCN~ etc.) is of the order of 10~ 2 , 
then the timescale for conversion is 

t ~ , 1 - NO ~ 1(T 12 s. (2) 

k x (10-V) 2 v ; 

Therefore, the chemistry is very rapid indeed. So long as this timescale is very much less 
than the dynamical timescale, then the details of the expansion and cooling of the gas are of 
minimal relevance; essentially the chemistry takes place very shortly after the sublimation 
of the ice mantles when the gas density is highest; the relative abundances of the product 
molecules are then "frozen in" to the flow. These qualitative conclusions are confirmed by 
the detailed results presented in Section 3. 

We consider a sphere of ice, instantaneously sublimated into the gas-phase, of initial 
radius r and density n = 10 23 molecules cm -3 which expands into a vacuum at some 
fraction, e, of the sound speed v s . The parameter e allows for deviations from completely 
free, spherically symmetric, expansion (e = 1) so that a value of e = would correspond to 
the situation of perfect trapping of the gas in a cavity. Thus, by mass conservation 

n o r o = m ' 3 = n ( r o + ev s t) 3 (3) 
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- = (-?-})'■ w 

no \r + ev s tj 

If r is assumed to be comparable to the typical thickness of an ice mantle on a dust grain, 
then r = 10~ 5 cm, v = 10 4 cms -1 and so the evolution of number density would in this 
case be given by 

n/n = 1/ (1 + 10 9 et) 3 (5) 

The initial temperature of the gas, To, is unknown as is the way that it varies with time. 
We are guided here by the observational results that indicate hot core temperatures of a 
few hundred Kelvin. The intial temperature may be higher, so we examine the temperature 
sensitivity up to 500 K. 

For an ideal gas expanding into a vacuum (Joule expansion) the temperature would 
be independent of time. For adiabatic expansion, TV 1 " 1 = constant. Thus, for spherical 
expansion of a gas with 7 = 4/3, T oc r _1 . In reality, the gas will cool as a result of work 
against intermolecular forces. In our model we consider limiting cases of 

T = To(^) 9 (6) 

where q — or 1. The various free parameters in our model are listed in Table [5] together 
with the range of values that we have investigated. For most of our model calculations we 
use the following values: T = 100, q — 0, a — 0, e = 1, with ko and tiq variable. We also 
generally assume that the third reactant is passive (A3 = off, see Table [5]) 

We calculate the time-dependences of the chemical models using a simple model, based 
around the GEAR integration package. Although the three-body reaction network and the 
extremely high densities are most unlike those applicable to models of molecular clouds, the 
principles are the same. As there are no variations in the chemical rate coefficients, the 
resultant sets of differential equations are not numerically stiff. 



3. Results 

The general characteristics of the results are illustrated in Tables [6] and [7] for the case 
of n = 10 21 cm" 3 , T = 500 K, k = 10~ 33 cm 6 s~\ T oc l/r (q = 1), e = 1, a = -2 and 
A 3 = off. 

As can be seen from Table 6, in this model the temperature falls from its initial value to 
200 K by about 10 -9 s. The parent molecules show a slight loss in relative abundance over this 
period and then remain constant. Conversely, we see in Table 7 that product molecules grow 
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rapidly in relative abundance during a period of less than 10~ 10 s and achieve steady-state 
within about 10~ 9 s. The steady-state abundances depend on the two-stage chemistry and 
on the relative abundances of the parent species. We emphasise that it is assumed that all 
rate coefficients are the same and have the same temperature dependences, so there is no 
selectivity in the chemistry other than that introduced by the network and the initial parent 
abundances. In this particular example, products Pi P5 attain abundances relative to 
H 2 of about 1CT 4 , products P 6 P15 reach about 10~ 5 , while products Pi 6 P 2 6 are 
in the range 10~ 2 — 10~ 3 . We note at this point that an abundance of 10 -5 relative to H 2 
may correspond to a column density in a typical hot core of about 10 15 cm -2 . We return to 
this point in Section 4. 

An important feature of the model whose results are given in Tables El and [7J is that 
the chemistry is dominated by reactions occurring at the very highest densities and earliest 
times. Therefore, the relative chemistry arising from these reactions is essentially unchanged 
during the later expansion and may be said to be "frozen-in" to the expanding gas. The 
dominance of the early time in the chemistry (at least, for this model) also means that the 
dependence of the gas temperature on the expansion time (or radius) is fairly unimportant; 
this is also the case for the dependence of the rate coefficient on temperature. For example, 
for the case reported in Tables |6] and [71 if all the temperature dependences are removed then 
the relative abundances are not significantly changed. In what follows, we shall ignore all 
temperature dependences. 

Clearly, the assumed initial temperature, T , may influence the model results. In fact, 
a change in To from 500 K (as in Table 2) to 100 K makes only rather slight changes in the 
relative abundances. In the model whose results are shown in Tables [6] and k oc T -2 . The 
most significant of the remaining parameters are the initial density and the rate coefficients. 
We have therefore explored the density/rate coefficient space, over a wide range. We plot 
in Figure 1 a summary of these results for the case where at least some of the product 
molecules have relative abundances of 10~ 5 or larger. There is a wide range of (n, k) space 
where chemistry meets this requirement for a trapping parameter of either 0.1 and 1.0. 

For the highest density case, some chemistry occurs on timescales of about 10~ 10 s even 
when the rate coefficient is as low as 10~ 41 cm 6 s _1 ; we find that products Pi6 - P 2 6 have 
relative abundances in this case of 10~ 6 — 10~ 7 . We interpret the rate coefficient for these 
small values by assuming that it has a form proportional to exp(— 0/T) (see equation [TJ), 
where /3 is a barrier height. Then 

(3 = Tln{k /k). (7) 

If ko = 10~ 31 cm 6 s _1 (a plausible value), then the barrier height for T = 100 K implied 
by a rate coefficient of 10 -41 cm 6 s _1 is 2300 K, a typical barrier height for bond-breaking 
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of an H-X bond (see Section 4).. For the lower densities, then the value of k required to 
produce a significant chemistry is of course much larger. For example, at an initial density 
of 10 20 cm -3 with zero temperature dependence, a value of k — 10~ 31 cm 6 s" 1 produces 
relative abundances of products Pi - P15 of 10~ 5 — 10~ 6 . A value of 10 -31 cm 6 s _1 implies a 
barrier height of less than 500 K, which may be implausible. At an even lower initial density 
of 10 17 cm~ 3 , then even the maximum plausible value of k, 10~ 29 cm 6 s _1 , fails to produce 
product molecule relative abundances larger than 10~ 6 . 

The effect of changing the gas expansion parameter, e, from its free expansion value of 
1 to a value of 0.1 that may represent a hindered expansion from a partially enclosed cavity, 
may also be considered. This reduction in e has the effect of maintaining the gas at the 
highest density for a longer time, so that the effect is allow the chemistry to proceed more 
quickly. For example, we compare the case of gas at a density of 10 22 cm -3 where all rate 
coefficients are 10~ 37 cm 6 s _1 . For this case, the product molecules range over four orders 
of magnitude in relative abundances. In all cases, the relative abundance of a particular 
product molecule is larger by one order of magnitude when the evaporation is restricted. 
This behaviour is also recovered for a case with lower density, 10 21 cm -3 in which some 
products are present with near-zero abundances. 

4. Discussion and Conclusions 

There are a number of general conclusions that may be drawn from this feasibility study: 

1. since the density falls rapidly in the evaporating gas, the chemistry is dominated by 
reactions occurring at the very earliest phases, typically within about 10~ 10 s, and 
remains "frozen-in" at later times; if reactions work - i.e. if activation barriers can 
be overcome - then the chemistry will be fast and efficient; 

2. as a consequence, any assumed temperature evolution in the gas is unimportant, and 
any assumed temperature-dependence in the adopted rate coefficients is also of little 
consequence; 

3. therefore, the dominating parameters are the initial gas number density, no, and the 
adopted three-body rate coefficients, k; 

4. the (n, k) parameter space explored includes a large region in which the complex chem- 
istry is rich enough to be observationally significant; i.e. where complex products have 
an abundance relative to water of about 10~ 5 , which could be relevant for hot core 
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chemistry. The region of (n, k) parameter space in which this occurs is indicated 
schematically in Fig.l; 

5. assuming that the form of the rate coefficients is 

k = k (T/300 K) Q exp (-/3/T) (8) 

where ft is the barrier height in Kelvins for the assumed reactions, then, if ko is assumed 
to have a canonical value of 10~ 30 cm 6 s _1 , the range in k of 10~ 30 — 10 -41 cm 6 s _1 
corresponds to a range in (3 of to 2530 K if the gas temperature is about 100 K; 
this range in /3 is plausible for simple bond-breaking reactions, if the bar rier height is 



roughly about 5% of the bond energy (IGlasstone. Laidler fc Eyringlll94ll ); 



6. an important conclusion is that there is a limit to chemical complexity imposed by 
the very short timescale available; for the range of parameters investigated here, this 
complexity corresponds to molecules containing about three or four heavy atoms (C, 
N, O); larger species are not expected from the mechanism described here; 

7. if the ionization level is even slightly enhanced from zero (by some suitable process, not 
considered here) then the rates an conversion efficiencies may be significantly higher; 
similarly, if radicals were to be introduced into the expanding gas ( a process not con 



sidered here) , then a chemistry of the kind originally envisaged by lAUen fc Robinson 



( 119771 ) may occur. 



The overall conclusion is this: we have demonstrated a proof of concept that a rich 
chemistry may occur in gas evaporating from a chemically-mixed ice in the co-desorption 
event, assuming that the event is sufficiently narrow in temperature for very high gas densities 
to be achieved in the evaporate. A similar conclusion may be made for other desorption 
events, such as the so-called 'volcano' event, but the composition of the evaporating gas and 
the products would be different. 

The results from this very preliminary work are of course higly tentative, but do seem to 
indicate, firstly, that there could be astronomical consequences from this idea, and, secondly, 
that further study of the idea (and of competing ideas involving three-body chemistry) 
probably cannot be made theoretically and would require laboratory investigation. 

We may crudely regard any chemistry that generates product molecules at about 10~ 5 
relative to H 2 as potentially important from the astronomical point of view. Hot cores 
are usually assumed to have about 1000 visual magnitudes of extincti on, corresponding 



appro ximately to a column density of hydrogen of about 10 cm (e.g., iMillar fc Hatchell 



19981 ). The available oxygen is mainly divided between CO (formed in the gas phase in 
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the early stage of the gravitational collapse of the star-forming core) and H2O (formed by 
hydrogenation of O-atoms at the surface of dust and retained there as the ice mantle). The 
column density of each species (i.e. CO and H 2 0) for hot cores in the Milky Way will be on 
the order of 10 20 cm -2 . Therefore, molecules with an abundance relative to water of 10~ 5 
may be expected to have column densities around 10 15 c m" 2 . The ty pical column density 
range of organic molecules in LMH is 10 15 - 10 16 cm 2 (jSnyde r 2006) ; for example a mino 
acetonitrile (NH 2 CH 2 CN) with a column density of 2.8 x 10 16 cm" 2 JBelloche et al.ll2008h. an d 
ethyl formate (C 2 H 5 OCHO) with a column density of 5.4 x 10 16 cm" 2 JBelloche et al.ll2009l) 
For an extended review of complex organics in interstellar clouds see iHerbst &; van Dishoeck 
(|2009h . 



Snyderl (120061 ) also emphasises that those molecules are hydrogen-poor. That seems 



likely to be the case for molecules arising in the processes considered here, as H-atoms 
will be ejected in the bond-breaking and bond-making processes of three-body chemistry. 
Indeed, the values of the rate coefficient that seem appropriate in this model do seem to 
imply bond-breaking of bonds of energies a few eV, i.e., corresponding to bonds involving 
H-atoms. We suggest that an experimental investigation be made to test the validity of the 
ideas expressed here. 
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Table 1: Initial gas-phase composition^, relative to H 2 (100) (jWhittetll2002l ) 
Species Abundance 



H 2 

CO 

CH 4 
H 2 CO 
CH 3 OH 

NH 3 
OCN- 



100 

15 

4 

3 

3 

1 

< 1 



( a )We do not include CO2 in the chemistry, as we regard that species as chemically inert in this context. 



Table 2: Postulated reactions between three active partners 
Number Reaction 



1 


CO 


+ 


CH 4 


+ 


H 2 


->• 


Pi 


2 


CO 


+ 


OCN- 


+ 


H 2 


-> 


P2 


3 


CO 


+ 


H 2 CO 


+ 


H 2 


-> 


P 3 


4 


CO 


+ 


CH3OH 


+ 


H 2 


-> 


P 4 


5 


CO 


+ 


NH 3 


+ 


H 2 


— > 


P5 


6 


CH 4 


+ 


OCN~ 


+ 


H 2 


— >• 


P 6 


7 


CH 4 


+ 


H 2 CO 


+ 


H 2 


— > 


Pt 


8 


CH 4 


+ 


CH3OH 


+ 


H 2 


-> 


P 8 


9 


CH 4 


+ 


NH 3 


+ 


H 2 


— y 


P 9 


10 


OCN- 


+ 


H 2 CO 


+ 


H 2 


-> 


P10 


11 


OCN- 


+ 


CH3OH 


+ 


H 2 


-> 


Pu 


12 


OCN- 


+ 


NH 3 


+ 


H 2 


-> 


P12 


13 


H 2 CO 


+ 


CH3OH 


+ 


H 2 


-> 


Pl3 


14 


H 2 CO 


+ 


NH 3 


+ 


H 2 


-> 


Pl4 


15 


CH3OH 


+ 


NH 3 


+ 


H 2 


-> 


Pl5 
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Fig. 1. — The region above the lines indicate the values of n and k for which some products 
have abundances larger than 1CT 5 relative to H 2 0, for a model in which T = 100 K, q = 0, 
and a — 0. When the expansion is partially hindered, e = 0.1, the chemistry is faster. 
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Table 3: Postulated reactions between two active and one passive partner: Stage I 



Number Reaction 



16 


CO 


+ 


H 2 


+ 


H 2 


-)■ 


Pie 


17 


CH 4 


+ 


H 2 


+ 


H 2 


->■ 


Pl7 


18 


OCN- 


+ 


H 2 


+ 


H 2 


-> 


Pl8 


19 


H 2 CO 


+ 


H 2 


+ 


H 2 


-> 


Pl9 


20 


CH 3 OH 


+ 


H 2 


+ 


H 2 


-> 


P20 


21 


NH 3 


+ 


H 2 


+ 


H 2 


-> 


P21 


22 


CO 


+ 


CH 4 


+ 


H 2 


-> 


P22 


23 


CO 


+ 


OCN- 


+ 


H 2 


— > 


P23 


24 


CO 


+ 


H 2 CO 


+ 


H 2 


-> 


?24 


25 


CO 


+ 


CH3OH 


+ 


H 2 


-> 


P 25 


26 


CO 


+ 


NH 3 


+ 


H 2 


->■ 


?26 
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Table 4: Postulated reactions between two active and one passive partner: Stage II 



Number Reaction 



27 


Pl6 


+ 


CH 4 


+ 


H 2 


->■ 


Pi 


28 


Pl6 


+ 


OCN- 


+ 


H 2 


-> 


P 2 


29 


Pl6 


+ 


H 2 CO 


+ 


H 2 


-> 


P 3 


30 


Pl6 


+ 


CH3OH 


+ 


H 2 


-> 


P 4 


31 


Pl6 


+ 


NH 3 


+ 


H 2 


-> 


P 5 


32 


Pl7 


+ 


OCN- 


+ 


H 2 


-> 


P 6 


33 


Pl7 


+ 


H 2 CO 


+ 


H 2 


— > 


P7 


34 


Pl7 


+ 


CH3OH 


+ 


H 2 


-> 


P 8 


35 


Pl7 


+ 


NH 3 


+ 


H 2 


-> 


P 9 


36 


Pl8 


+ 


H 2 CO 


+ 


H 2 


-> 


P10 


37 


Pl8 


+ 


f A T T / \ T T 

CH3OH 


+ 


H 2 


— > 


Pll 


38 


Pl8 


+ 


NH 3 


+ 


H 2 


-> 


Pl2 


39 


Pl9 


+ 


CH3OH 


+ 


H 2 


-> 


Pl3 


40 


Pl9 


+ 


NH 3 


+ 


H 2 


-> 


Pl4 


41 


?20 


+ 


NH 3 


+ 


H 2 


— > 


Pl5 


42 


P21 


+ 


CO 


+ 


H 2 


-> 


P 5 


43 


P21 


+ 


CH 4 


+ 


H 2 


-> 


P 9 


44 


P21 


+ 


OCN- 


+ 


H 2 


-> 


Pl2 


45 


P21 


+ 


H 2 CO 


+ 


H 2 


-> 


Pl4 


46 


P21 


+ 


CH3OH 


+ 


H 2 


-> 


Pl5 


47 


P22 


+ 


H 2 


+ 


H 2 


-> 


Pi 


48 


P23 


+ 


H 2 


+ 


H 2 


-> 


P 2 


49 


?24 


+ 


H 2 


+ 


H 2 


-> 


P 3 


50 


P 25 


+ 


H 2 


+ 


H 2 


-> 


P 4 


51 


?26 


+ 


H 2 


+ 


H 2 


-> 


P 5 
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Table 5: Free parameters in the model 
Parameter Description Value 



k Rate coefficient 10~ 29 - 10" 41 cm 6 s" 1 

a Temperature dependence of rates (k oc T a ) or —2 

n Initial number density 10 17 — 10 23 cm -3 

e Trapping parameter 0.1 or 1.0 

T Initial temperature 100 — 500 K 

q Radial dependence of temperature (T oc r~ q ) or 1 

A 3 Flag to set reactions involving 3 active partners on or off 



Table 6: Results from the model described in the text. Results are shown for the reactant 
species and are given as the log of the species abundance relative to H 2 0. 



Time (sec) 


Density (cm 3 ) 


Temp (K) 


CO 


CH 4 , HCN 


H 2 CO, CH 3 OH 


NH 3 


7.98xl0- n 


7.89 xlO 20 


462.1 


-0.86 


-1.43 


-1.55 


-2.03 


1.66xl0- 10 


6.26xl0 20 


427.7 


-0.88 


-1.45 


-1.58 


-2.05 


2.59xl0~ 10 


4.95 xlO 20 


395.4 


-0.89 


-1.46 


-1.59 


-2.07 


3.59xl0~ 10 


3.94xl0 20 


366.6 


-0.90 


-1.47 


-1.60 


-2.08 


4.68xl0- 10 


3.09xl0 20 


337.9 


-0.90 


-1.48 


-1.60 


-2.08 


5.85xl0- 10 


2.51 xlO 20 


315.4 


-0.91 


-1.48 


-1.61 


-2.09 


7.11xl0- 10 


1.99xl0 20 


292.1 


-0.91 


-1.49 


-1.61 


-2.09 


8.48xl0~ 10 


1.55xl0 20 


268.7 


-0.91 


-1.49 


-1.61 


-2.09 


9.95xl0~ 10 


1.24xl0 20 


249.7 


-0.91 


-1.49 


-1.61 


-2.09 


1.15xl0- 9 


9.40xl0 19 


227.3 


-0.91 


-1.49 


-1.61 


-2.09 


1.33xl0- 9 


7.78xl0 19 


213.4 


-0.91 


-1.49 


-1.61 


-2.09 


1.51xl0- 9 


6.05xl0 19 


196.3 


-0.91 


-1.49 


-1.61 


-2.09 
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Table 7: Results from the model described in the text. Results are shown for the product 
species and are given as the log of the species abundance relative to H 2 0. 



Time (sec) 


Pi, P 2 


P 3 , P 4 


?5 


P 6 


P 7 , P 8 , P 10 , Pn 


P 9 , P12 


P13 


Pu, P15 


7.98xl0 -11 


-4.66 


-4.78 


-5.08 


-5.53 


-5.66 


-5.83 


-5.78 


-5.96 


1.66xl(r 10 


-4.34 


-4.36 


-4.66 


-5.11 


-5.23 


-5.41 


-5.36 


-5.53 


2.59xl(r 10 


-4.04 


-4.17 


-4.47 


-4.91 


-5.04 


-5.22 


-5.16 


-5.34 


3.59XKT 10 


-3.93 


-4.06 


-4.36 


-4.81 


-4.93 


-5.11 


-5.06 


-5.23 


4.68XKT 10 


-3.87 


-3.99 


-4.30 


-4.74 


-4.87 


-5.05 


-4.99 


-5.17 


5.85XKT 10 


-3.83 


-3.95 


-4.26 


-4.70 


-4.83 


-5.01 


-4.95 


-5.13 


7.11x10 iU 


-3.80 


-3.93 


-4.23 


-4.68 


-4.80 


-4.98 


-4.92 


-5.10 


8.48xl(T 10 


-3.78 


-3.91 


-4.21 


-4.66 


-4.78 


-4.96 


-4.91 


-5.09 


9.95xl0- 10 


-3.77 


-3.90 


-4.20 


-4.65 


-4.77 


-4.95 


-4.90 


-5.07 


1.15xl0~ 9 


-3.76 


-3.89 


-4.19 


-4.64 


-4.76 


-4.94 


-4.89 


-5.07 


Time (sec) 


Pie 


Pl7) Pl8 


Pl9) ?20 


P21 


P22, P23 


P24, P25 


P26 




7.98XKT 11 


-2.05 


-2.62 


-2.74 


-3.22 


-3.47 


-3.59 


-4.07 




1.66xl0- 10 


-1.84 


-2.41 


-2.53 


-3.02 


-3.27 


-3.40 


-3.87 




2.59xl0~ 10 


-1.74 


-2.31 


-2.44 


-2.92 


-3.19 


-3.31 


-3.79 




3.59xl0- 10 


-1.69 


-2.26 


-2.38 


-2.87 


-3.14 


-3.27 


-3.75 




4.68xl0- 10 


-1.66 


-2.23 


-2.35 


-2.84 


-3.12 


-3.24 


-3.72 




5.85xl(T 10 


-1.64 


-2.21 


-2.33 


-2.82 


-3.10 


-3.22 


-3.70 




7.11xl0- 10 


-1.62 


-2.20 


-2.32 


-2.81 


-3.09 


-3.21 


-3.69 




8.48xl0~ 10 


-1.62 


-2.19 


-2.31 


-2.80 


-3.08 


-3.21 


-3.68 




9.95xl0~ 10 


-1.61 


-2.18 


-2.31 


-2.79 


-3.08 


-3.20 


-3.68 




1.15xl0~ 9 


-1.61 


-2.18 


-2.30 


-2.79 


-3.07 


-3.20 


-3.68 





